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ABSTRACT 


The theory and data reduction procedure for constant temperature 
hot wire anemometry are presented. The procedure is valid for all 
Mach and Prandtl numbers, but limited to Reynolds numbers based on 
wire diameter between 0.1 and 300. The fluids are limited to gases 
which approximate Ideal gas behavior. Losses due to radiation, free 
convection and conduction are included. 
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NOMENCLATURE 


Most of the nomenclature Is Ms ted; those excepted are the 
numerous factors In the correlation equations which are obvious In the 
context. 

a non-dimensional temperature 

A a(£) - a* 

C N Nusselt number end loss correction factor 

C R recovery ratio end loss correction factor 

d wire diameter (In.) 

2 

g gravitational acceleration, 32.17 ft/sec 

Gr Grasshof number 

p 

h convective heat transfer (BTU/hr-ft -°R) 

I current (amps) 

k thermal conductivity (BTU/hr-ft-°R) 

K unit length resistance (ohms/in.) 

Kn Knudsen number 

X, effective wire length (in.) 

M Mach number 

Nu Nusselt number 

2 

p pressure (lbf/in ) 

Pr Prandtl number 

q heat transfer rate (BTU/hr) 

q m measured hot wire power dissipation (watts/sec) 

r wire resistance (ohms) 

R gas constant (ft-lbf/lbm °R) 

S end loss correction parameter 

111 


T temperature (°R) 

U velocity (ft/sec) 

x length coordinate (in.) 

z non-dimensional length coordinate 

temperature coefficient of resistance, (°R ) 

o 

a 2 second temperature coefficient of resistance (°R ) 

6 Isothermal expansivity coefficient (°R _1 ) 

Y specific heat ratio 

e emisslvity 

C temperature difference, T(x) - 7\ 

n temperature recovery ratio, T e /T s 

n normalized recovery ratio 

n c recovery ratio at continuum limit 

rif recovery ratio at free molecule limit 

0 flow direction angle 

p viscosity (Ibm/ft-hr) 

£ non-dimensional length coordinate 

p density (lbm/ft 3 ) 

cx Stefan-Boltzmann constant 

to recovery ratio end loss correction parameter 

Subscripts 


( 

> a 

air 

( 

>c 

conduction 

( 

>d 

duct or tunnel walls 

( 

>e 

equilibrium 

{ 

>f 

fluid 

( 

) m 

measured 



t 


( ) r radiation 

( ) s wire support 

( ) t total or stagnation 

( ) w1 wire 

{ ) 0 reference values except on a 

( ) local free stream if on M, Re, or e; infinite wire if on Nu 

1 '00 

( ) average 


v 
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1.0 INTRODUCTION 

The theory and data reduction procedure for hot wire anemometry are 
presented in this paper. The basis of hot wire anemometry is that the 
heat transfer rate is a function of the flow properties of the flow to 
be measured. By heating the wire above the equilibrium temperature and 
measuring the current fluctuations to maintain the wire at a constant 
temperature, the heat loss and hence the flow speed is determined via 
a series of Mach, Reynolds, and Nusselt (heat transfer coefficient) 
correlation curves. 

The hot wire loses heat by all three heat transfer modes: convection, 
conduction, and radiation. The effect of each one is considered, since 
the correlation curves are valid only for the forced convective heat loss. 
Due to the finite length of the wire, the conduction end losses are 
considerable, which affects both the heat transfer (convective) and the 
cold wire equilibrium temperature, which again affects the heat transfer. 
An analysis of all these "end" losses, free convective and radiative 
losses is given. 

A complete data reduction procedure is given. The computer program 
and its data input requirements are given. The appendix also presents 
the non-linear effects of temperature on the wire resistance and a 
limited uncertainly analysis of the Mach number due to uncertainties in 
the wire diameter. 

The technique for obtaining the data is not presented here. The 
program listed is quite general, being valid for all types of gases at 
all Mach numbers for Reynolds numbers based on wire diameter between 
0.1 and 300. 
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2.0 HOT MIRE ANEMOMETRY THEORY 

The hot wire probe is a device for measuring the local mass flux of 
a flowing fluid. The probe consists of a very fine wire electrically 
heated and immersed In a fluid. By measuring the heat transfer rate off 
the wire, the mass flow can be deduced. The heat transfer rate is easily 
obtained by measuring the electrical resistance and the voltage drop of 
the hot wire. 

Two principles are the basis of hot wire anemometry: 1) the 
electrical resistivity of the wire is temperature dependent, preferably 
linearly, and 2) the heat transfer rate is a function of the flow 
properties of the flow. The first principle is necessary so that the 
temperature of the wire needed for the heat transfer determination can 
be accurately measured. The need for the second principle is obvious. 

The hot wire has two modes of operation: constant current and con- 
stant temperature. In the constant current mode, the electrical current 
through the wire is held constant by an electronic circuit. A variation in 
the mass flow changes the heat transfer, temperature and, thus, resistance 
of the wire. The instantaneous resistance is measured to obtain the heat 
transfer rate and the mass flow. In the second mode the current is modu- 
lated with the varying flow to maintain it at constant temperature and 
resistance. The instantaneous current through the wire is measured to 
obtain the heat loss. The constant temperature mode has several advan- 
tages over the constant current mode, among which are simpler and more 
straightforward data acquisition and reduction procedures. These and 
other reasons have made the constant temperature method more useful and 
will be used in this report. 

The theoretical analysis of heat transfer from circular cylinders 
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will not be given here. As yet no complete analysis with all the compli- 
cations exists. Partial theories have been developed for the free molec- 
ular, high Reynolds number, and low Reynolds number flow regimes. For 
the transition zones empirical data only are available. Extensive empirical 
data have been published on the heat transfer rates in terms of Nusselts 
number of infinitely long circular cylinders as a function of the flow 
properties, namely Mach and Reynolds numbers. These data are reduced into 
a set of "universal" correlation curves, Fig. 3. The universal correlation 
curves allow the direct determination of the Mach and Reynolds (based on 
wire diameter) from the measured heat transfer corrected for the effects 
of finite wire length. 

The finite wire differs from the infinitely long wire in that heat 
is conducted to the supports of the finite wire in addition to the con- 
vected and radiated heat loss. This conduction end loss causes a nonuniform 
temperature distribution. The infinite wire loses heat only by convection 
and radiation and therefore has a uniform temperature distribution. The 
end loss of the finite wire also affects the adiabatic recovery temperature 
of the wire which also affects the heat transfer rate of the finite wire 
as compared to the infinite wire. 

The end loss effects of the finite wire must, therefore, be carefully 
taken into account before the correlation curves of Fig. 3 for infinitely 
long cylinders are applicable. Since the end effects are so important 
(the end conduction heat loss can be as great as 50% of the total heat loss 
for wire aspect ratios of 100), the theoretical analyses of the end losses 
are given in detail below. 

2.1 Derivation 


The system to be analyzed is shown in Fig. la and b. Figure la 
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shows the complete wire probe and lb, the wire and coordinate system. The 
steady state case will be considered only. 

The wire segment of length Ax loses heat by four mechanisms, which are 
forced convection 

% - h » d Ax'T(k)-T e ) 

where h = coefficient of convective heat transfer 
T(x) = temperature of wire 

T fi = equilibrium (adiabatic wall) temperature of wire, 

free convection 


$f.c. ■ h f.c. » <1 *<T<m-V 

where h* „ = average free convective heat transfer coefficient 

~ 1.16 (Gr Pr) 1/10 (Ref. 21) 

for 10“ 5 < Gr < 10 3 

p 2 ge(T(k)-TJd 3 

and Gr = Grashoff *s number = * — ~ — » 

Vf 


radiation 


q r * a e tt d Ax(T(k) 4 -T d 4 ) 

where T d - duct temperature* 


and 

conduction 




d 2 T(x) 

dx 


4 vV haat losses are balanced by the ohmic heating of the wire 


^generated * 1 2 Ax 

where I = electrical current through the wire 
and 

r(x) = resistance per unit length of the wire at 
temperature T(x) 

* r 0 [l +o,(T(x)-T 0 ) + o 2 (T(x)-T 0 ) 2 ] 

T q = reference temperature 

r Q = reference resistance at T Q 

= first temperature coefficient of resistivity 

a 2 = second temperature coefficient. The Influence 
of otg is negligible for small overheat ratios 
C(r(x)-r e )/r e L as Is shown in Appendix B. 

Defining ?(x) = T(x) - T and assuming that T d ^ T fi and 
|[T(x)-T d )]/T d ]| < 1 or 

T(x) 4 - T d 4 * 4 T(x) 3 [T(x)-T e ] 


and simplifying yields the energy equation 

k 


h E(x) + 4 a e T(x) 3 c(x) + 1.16 — ^ (Gr Pr)^ 10 ?(x) 
kjrt d d 2 c(x) I 2 r c 


dx 


Td“ (1 + “1 ?(x)) ■ 


The boundary conditions are 
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Introducing the nondimensional variables 
K * and a(£) » ^ c(x) 

simplifies the differential equation to 


[M + UpTfoK ,, , )6(Gr pr) i/io 

T • K f 

= k wi d 2 d 2 a(£) + a l 1 2 r e 

"^7 7 dr 


0(i 



] a(5) 


The first and last terms within the brackets (forced convection and ohmic 
heating) are of 0(1). The other two terms within the brackets (radiation 
and free convection) are of smaller order as shown and hence are negligible 
Typical values for air flows at M = 2.5 with overheat ratios of 1 and 
d - 0.0004" are T(£) * 1000° R, T e * 500° R, k f * 0.02 BTU/hr-ft-°R, 

C p = 0.25 BTU/lbm-°R, y f = 0.05 Ibm/hr-ft, p f = 0.006 lbm/ft 3 , 
g = 32 ft/sec 2 and 0 = 0,002/°R. The radiation term 


4 q_e T?_(£)_d = 10 -2 
k f 

and the free convection term 

1 ,16(Gr Pr) 1/10 = 0.2. 


Therefore radiation and free convection are negligible provided that the 
overheat ratios are less than unity and wire diameters do not exceed 
0.0004". If they are larger the free convection term must be considered. 
A theoretical discussion of free and forced convection of fine hot wires 
is given by Woods (22). 

The governing differential equation reduces to 
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2 

^a( 5 ) (, ta(5) ) 

with the boundary conditions 
i||2l = 0 and 


wi d 

"T 7 


£ a U ) 

A'C 


a(±l) = o,(T s -T e ) 


2.2 Nusselt Number Correction 

The relation between the Nusselt humber (hd/kf) 
infinitely long wire and a finite wire must be determined, 
the infinite wire bj ( )* 

★ . * 

a (5) = a = constant, i.e., no end losses. 


Therefore the equation becomes 

* 


h = 


a I 2 r 


it d 


e 1 + a 
IT 


The average convective coefficient of a finite wire is 


F « 


a I 2 r 


It d 


e 1 + a 


where a - = average non-dimensional temperature of 
wire 

and also 

„ r -r 

a = V T w"V = (overheat ratio) 


Nu 
00 



h 


F 


, * — 
1 + a a 

1 + F a 


for an 
Denoting 


the finite 


Therefore 
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which gives the Nusselt number of an infinite wire in terms of the measured 

i 'ff _ ft 

(average) Nusselt number and a/a . To determine a/a , the differential 
equation is solved for a(£) and averaged over the wire length. 

On substituting for h and a-jl r e /rrd, in terms of F, a", and a , 
and rearranging > the differential equation becomes 


Fd 

TF 


1 + a 


a 

T 

a 


(a(5) - a ) * 


*wi 


d 2 

7 


^4 

df 


r\ 

hJ.. 


The measured Nusselt number Nu m = Fd/k^ and letting 
S * d/£ [(l+a/Nu m )(k wi /k f )] 1/2 obtains 

1 .^ 1 ^.^) +a * = 0 
a dr 


which becomes with one more change of variables 
z * C I* (a/a*) 1 ^ 2 and A(z) = a(£) - a* 

_ a(z) * 0 


with the boundary conditions 

= o and A[± (a/a*) 1/2 (1/S)] = (T s ~T e ) - a* 


The solution is 


A(z) 


VW - a* 
cosht(VS)(T/a*) ,/2 ] 


cosh (z) 


or 
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a( 5 ) = a* + cos ht i^^V^ A [o ,( T .-T e ) - a*] 
coshCfl/S) (a/a ) 1/z ] 


Integrating with respect to 5 and averaging yields 
[1 - (T -T )] 

4 = ! a ! r 1 /T/ a *\ 1/2 

a 


i <a/a *) 1/2 


tanh [ (a/a *) 1/f2 ] 


Kovasznay (1) assumed that the wire support temperature T s equals the 
cold wire equilibrium temperature T D and also that is unaffected 
by the end losses. With this assumption, the end loss correction factor 
for the Nusselt number is 

a tanh[ 3 - (a/a*) 1 / 2 ] 

7 = "T<iyIW " 


The graphical solution of this relation is shown in Fig. 2a. 

Dewey (10) included the effects of different wire support temperatures. 
These effects are also shown in Fig. 2a and are seen to be secondary so they 
can be neglected. 

The equilibrium (adiabatic wall) temperature of finite wire as measured 
differs from that of an infinite wire due to the end heat conduction losses. 
This difference, Ignored by Kovasznay (1), has been considered by Dewey (10). 
The equilibrium temperature ratio or recovery ratio end loss correction as 
a function of S is shown in Fig. 2b. 

The recovery temperature correction requires knowledge of the flow 
properties of the flow which are as yet unknown. Therefore in this paper 
the two corrections are considered separately. The Nusselt number 
corrections are given by Kovasznay in Fig. 2a and the recovery temperature 
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corrections by Dewey In Fig. 2b. Applying the Nusselt number corrections 
with the measured recovery temperature, a first approximation of the flow 
properties, Mach and Reynolds numbers, is obtained from the heat transfer 
correlation curves described below. These then can be used to correct for 
the recovery temperature end loss effect. The corrected recovery temper- 
ature is then used in place of the measured recovery temperature and the 
entire procedure is repeated. The procedure for correcting the recovery 
temperature will be described after the discussion of the heat transfer 
correlation curves. 

2.3 Nusselt -Mach-Reynolds Numbers Correlation 

Theoretical analyses for predicting the heat transfer to infin- 
itely long cylinders in the continuum flow and free molecule flow are 
available, but the transition between the two regimes must be obtained 
empirically. Since extensive empirical dat* have been published, no 
theoretical analyses will be presented here. The experimental data have 
been summarized and correlated by Dewey (10) in terms of Nusselt number 
of infinitely long cylinders as functions of the local free stream Mach 
and Reynolds numbers. These correlations are presented graphically and 
analytically. 

The best correlations were found to exist if all the fluid properties 
were evaluated at the total free stream temperature, Laufer and 
McClellan (9). The correlation curves are shown in Fig. 3. The correla- 
tion equations are given in two parts, as shown, 

Nu„(Re„,MJ = NiijRe.,.) *(Re„,Mj 

where Nu (Re^, 00 ) is the Nusselt number of an infinite cylinder 
at M^ = “. (Subscript «> on Nu refers to infinite cylinder and on 




Re and M refers to the local free stream conditions. Re 
wire diameter.) 

NujRe^) = Re oo n {0.14 + N ] + N 2 ) 


n = 1 


Re 


0.6713 


5.142 + 2 Re. 


~OTfT 


0.2302 Re 


0.7114 


"OT7TT4 


15.44 + Re, 


No = ( 


0.01569 


0.3077 + Re 


0.7378 


)( 


15 

15 + Re. 


and 4(Re ,M ) is the deviation of the Nu at from the 

’ ’ txi OO* OO 


= 1 + <f>i * 4>2 * 4>3 


♦l 

<P 2 

4*3 


0.06309 


+ 0.5701 [( 


M. 


1 .222 


1 + M 1 ^ 


) 


1.569 


1] 


= 1.834 - 1.634 ( 


Re 


1.109 


2.765 + Re 
Re. 


T7T09 


, + [0 .300 - ^ ( )] 


These correlation equations are valid only in the range 
0.4 < M < “ and 0.01 < Re < 1000. 

•— DO — OO — 

For boundary layer investigations, correlations for M m < 0.4 
hence the following correlation equations were derived for 

0 < M m < 0.4 
— 00 


11 

based on 


Nu at = °°. 


are required. 
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The M m = 0 correlation, as obtained from the Oseen solution of Cole 
and Poshko (16) and the data of Hi 1 pert (17) is 

Nu oo( Re oo»°) = °* 19 + °- 66 Re oo°’ 45 

In the Mach range 0 < M^ < 0.4, a weighted logarithmic average 
between the M^ = 0 correlation above and the Dewey correlation at 
M « 0.4 is used. 

CO 

(F £og[Nu (Re wf 0)] + (1-F) *og[Nu„(Re„,0.4)]} 

NU - 10 

00 

where 

F = [0.618 + (2.5 MJ 0,9 ]" 1 - 0.618 
The extended correlation is valid only for the range 


0 £ M m < 0.4 and 0.1 < Re w < 300 


Some of the original data, Dewey's correlation, and the new correlation 
are shown in Fig. 3. 

Other correlations have been used, specifically the Yanta correlation (18) 
as used by Sfeir (3). Yanta modified the Dewey equations to extend M^ down 
to 0.1. However, the correlation equations as reported by Yanta and Sfeir 
produced no correlation at all for M^^O.3. This difficulty and the desire 
to extend the Nu m - M_ - Re m correlations to all Mach numbers required the 
addition to the Dewey equations for < 0.4. 

Measuring the local free stream static pressure, in addition to the hot 
wire measurements, the Reynolds number is obtained as a unique function of 
the Mach number 


Re. 





It can be seen then that the Mach number must be determined iteratively 
from the correlation equations. 
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2.4 Recovery Temperature Correction 

To obtain the first approximation of the Mach and Reynolds numbers, 

the measured recovery temperature, T . for the finite wire, was used in 

cm 

the procedure of the previous section. The correlation data are, however, 
valid only for Infinite wires, hence the recovery ratio or temperature as 
measured must also be corrected for the end loss effects. 

The procedure for correcting for the recovery ratio end loss as outlined 
by Dewey (10) Is as follows: 

The end loss recovery correction factor is 



- 2 *- - [1 - 4 r 1 HO -*)' 1 
% em 

- Cl - «( 5 s - HO -*)' 1 


where 

* 

T om = cold wire recovery temperature of the Infinite wire 
cm 

to = a parameter that depends primarily on S and is 
given in Fig. 2b, 
and 


n = T/T t , the temperature ratio or recovery ratio 
depending on T. 

Since it is not feasible to experimentally determine the wire support tem- 
perature, a good approximation is 

2 


• T 4 


1 + p r 1/2 • Xjjl M 

C « 
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A normalized recovery ratio n of a conducting cylinder Is defined 
in terms of the free molecular and high Reynolds number continuum limits 
recovery ratio, and n c . respectively, 

n - (n-n c )/(n f -n c ) 


The continuum limit and free molecular recovery ratios have been correlated 
by Morkovin (7) as 


n = 1 - 0.050 ( 


M. 


3.4 


1.175 + M 


XT 


00 


and 


rif - n c c 0.2167 ( 


M 


2.80 


0.8521 + M 


XW 


The normalized recovery ratio is empirically correlated in terms of the 
local free stream Knudsen number in Fig. 4. The equation of the universal 
curve shown in the figure Is 

Kn, 1 - 193 

11 * 0.4930 + Kn ' ,193 ' 

00 

and 

Kn^ = t (tty/2) (M^/ReJ] 1/2 

'ft 'it 

Sufficient information exists to determine (* T em /T t ) given 

the and Re^ from the first Nu correction and the Nu - M - Re 

correlation. The final Nu . M I, Re^ and T * are determined by 

00' 00 00 gm * 

repeating the Nu correction with the computed (or T em ) instead 
of the measured n m (or Tj. Usually less than five iterations are 
required to obtain 0.01% accuracy in n m . Sfeir (3) claims the same 
accuracy with only three iterations. 
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2.5 Local Flow Direction 

The Mach number obtained In the preceding section Is the actual 
Mach number only If the hot wire is perpendicular to the flow. If the wire 
Is not oriented perpendicularly to the flow then the Mach number obtained 
is only the component of the flow normal to the wire. By obtaining the 
Mach number components at two different wire orientations, the resultant 
Mach number and direction are obtained as follows: 

The resultant local flow direction relative to some reference line 
(e.g., a wind tunnel centerline as in Fig. 5) is given by 



Mj sin 
Mj cos e 2 


M 2 sin Of 
M 2 cos 


) 


and the resultant Mach number is 
M w « M^/cosf | - 0 1 + u) 

= m 2 /cos( | - e 2 + e) 

With the resultant local free stream Mach number, the remaining flow 
properties may be obtained 

t. * y 0 + *r w. 2 ) 

Poo = P «/ RT ® 
and 

This completes the basic hot wire anemometry theory. A more complete 
analysis required for the data reduction procedure would have to include 
an uncertainty analysis. No attempt was made to determine the uncertainty 
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except that due to the wire diameter. The wire diameters were very diffi- 
cult to measure accurately since they are on the order of 0.0004 to 0.0001 
inches. The manufacturer’s nominal values are used and are ±3% uncertain. 
The derivation of the uncertainty due to the wire diameter Is given in the 
appendix. 
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3.0 COMPUTATIONAL PROCEDURE 

In this section, the step by step computation procedure for reducing 
the hot wire data is outlined. The computer program (the listings are in 
the appendix) follows this procedure exactly. The procedure for obtaining 
the data will net be discussed. 

1. Calibrate the wire to obtain the first temperature resistivity 
coefficient, , 

(r 2 -r,) 

1 = FfTTpTjT 


Handbook values are used for the second temperature coefficient of 
resistivity, a 2 . The temperature of the wire is 


T = T, + — 

la. 


r-r. 


a. 


r-r 


i ) 2 


a. 


2. Determine effective wire length, £, from the manufacturer's unit 
length resistance at a reference temperature T^ and the actual 
measured resistance at T^ 





3. Also determine the wire diameter or use the manufacturer's nominal 
specification (see step 25 for the wire calibration). 

4. Start the run and determine the cold wire measured equilibrium temper- 
ature, the heated wire temperature, and the total temperature, T^. 
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em 


m 


T 1 + a 


1 r em~ r l 


a 2 , r enf r l ,2 


T + JL r *~ rem °2 

ern a l r em «, 


l C / CHI 

1 r l 

a 2 ( 

~ *TT ' r 


w r em_ j2 
em 


5. Determine the viscosity and conductivity of air at the total temperature. 
Sutherland's viscosity equation (Ref. 12) 

y t x 10 7 * 7.323 T 3/2 /(T + 198.72) ( ) 


and conductivity 


. = 1.1408 x IQ" 3 T 3/2 ( BTIJ \ 

a 7 + (441 .72/10^ * 6/ ^) lir-tt-*!* 


6. Determine the thermal conductivity of the wire as a function of temper- 
ature (Ref. 14) 


k wi ■ w ( t- t i>< w v; :% <,2 > 


and compute k w1 at 1 = 1^. 

7. Measure the power dissipation in the wire 

..2 


%i 1 r wm 


wm 


where v = voltage drop across the hot wire and Is given In terms of 
the bridge voltage, V, by 


wm 


^bridge + r lead + r wm^ 


v 


V 
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8. Assume an equilibrium recovery factor or use the measured value 
(rijp *- T em /T t ), and compute 


em n m ' t 


9. Compute the air conductivity from step 5 with T = T, 




10. Compute the measured Nusselt number 


Nu m “ 5R 


4.097 q 


wm em' "am 


11, Calculate the end loss parameter, S, 


c = ± ( kwi J_ !wm \ }/Z 

A Nu m r em 


and the overheat ratio 


r -r 
'wm em 


12. Read off the value for 


= f(S) 


from Fig. 2a or use the cubic polynomial 


= 1 - 1.100 S + 0.0999 - :.250 S k 


hence a = a w /(a w /a w ) 
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arid the Nusselt number correction factor for the end loss 


a w , 1 + a w v 
C N ■ -T ( — r- ) 

\ 1 + a w 


Thus the corrected Nusselt number for an Infinitely long wire Is 


Nu <* = C N Nu m 


13. Assume M and with the measured p , compute Re 




Re = 68.067 

V RT t> 


jy? MJyO + ^ M ^ 2 )] 1 ^ 2 


14. With these values for and Re^, determine Nu^ from the corre- 
lation Nu - Ma - Re correlation curves as given In Fig. 2 and 
Section 2.3. 


Nu m = ftRe^.Mj 


15. Compare ■‘■‘ils Nu^ with the Nu^ obtained In Step 12. Readjust 
In Step 13 and recompute until 

Nujstep 12) = Nujstep 14). 

16, Knowing the and Re^, compute the Knudsen number 

Kn » = ( V ) Kef 


the continuum recovery ratio 


M 


3.5 


(1.175 + 


» 


1.0 - 0.05 
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the free-molecule flow recovery ratio 


1 n c + 0.2167 ( 


M, 


2.80 


0.8521 + M 


TM 


and the universal recovery factor 
1 .193 

TT5d * 


Kn 


0.4930 + Kn 


17. The recovery factor for the infinite wire is then 


n s n (n f -n c ) + n c 


18. Compute the adiabatic wire support temperatures from M 


T s • T t (l + (Pr) 1/2 X£-mJ)/(1 + 


and the support temperature ratio 


'Is “ V T t 


19. Read the ratio w/1-id as a function of S from Fig. 2b or compute 
from 

w = 0.700 S + 1 .125 S 2 + 4.250 S 3 


and solve for w 


(0 « 




m 


20. The end loss recovery factor correction ratio is 


[1 - wfrig/r^)] 




OJ 
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21 . The "measured" wire recovery ratio Is now computed 

V * ^ C R 

22. Compare with the measured or assumed value of of Step 8. If the 
two differ, then go back to Step 8 and repeat the entire procedure 
with the computed ''measured' 1 wire recovery ratio n^'. Usually no 
more than five iterations are required to obtain ±0.01% precision 

in M w . 

23. The normal component has now been determined for one wire orienta 
tion. Repeating the entire procedure obtains another normal 
component at a different wire orientation at the same point in the 
flow field. With the given wire orientation angles 6 1 and e 2 and 
the computed normal components, the resulting H* direction and 
magnitude are 


0 = 


^ Mj sin e 2 - M 2 sin e-j 
tan ^ M-j cos 0 2 - M 2 cos 0<j ^ 


and 




= M-j/cos( \ - 8 1 + 6) 


24. The remaining flow properties are now computed from the resultant 

t» - y 0 + ^ n. 2 ) 

144 p„ 


and 


3 co * TV 




1/2 


23 


25. If the wire diameter used in Step 3 is of doubtful accuracy, the 

wire diameter is calibrated by measuring a very accurately known 

M w and adjusting the wire diameter until the computed equals 

the known M . 

00 

This completes the procedure for reducing the hot wire data at one 
particular point in the flow field. 

A final word of warning is necessary. Tha whole point of using the 
Nu - Ma - Re correlation curves is to eliminate the need for calibration 
curves. This avoids the danger of calibrating bad data into "good" data. 
Calibrating the wire diameter as in Step 25 reintroduces that danger, and 
any adjustments of wire diameter greater than the manufacturer's tolerances 
should be considered suspect. In other words, the source of any discrepancy 
lies elsewhere. 

The following section outlines the input data required for the data 
reduction program listed in the appendix. 
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4.0 DATA INPUT 

The data input required for the data reduction program can be divided 
into three types: program control data, general tunnel and wire conditions, 

and the specific data output of the hot wire anemometer. All data are 
initialized either at zero or a default value, as shown below. The data 
read-in is controlled by Subroutine DATA (IEOF), shown in Appendix A. 

Each group of data within each of the three types is identified by 
a preceding card containing an integer 1 through 20 in the first two 
columns, right justified. Thus there are 20 data groups. The only excep- 
tion to this rule is the data group labeled by the integer 1. This card 
must initially contain another integer in Columns 3 to 5, right justified. 
The other 19 identifiers contain one integer only; Columns 3 to 5 must be 
blank, otherwise the program stops. 

Decimal data following the identifier card must be formatted in six 
words per card of 10 columns per word and must be right justified if in 
exponential form. Again the exception to this rule is the data labeled 
by Integer 1. Here the format is 8 columns for the first word, and 10 
columns for each of the remaining 5 words. The reason for this exception 
is required by the peculiarities of the remote teletype terminal used to 
record the data. There is another difference in the Integer 1 data. The 
card immediately following the label card must either be blank or contain 
alpha-numeric characters in Columns 1 - 72. The same characters are 
printed in the output as a title. Omitting this card will cause a program 
termination. 

Integer data are only 2 words per card and must be right justified in 
Columns 1 - 2, and 3 to 5. There are only two groups requiring integer 
data, namely, the data labeled by integers 14 and 15. The data sets 


J 
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labeled by integers 19 and 20 are blank and have been provided for future 
expansion of the program. 

The three types of data required are described below, along with 
their default or initial values. 


Program Control Data 

This type of data is used to control the accuracy required, 
maximum number of iterations allowed, number of data sets to 
be reduced, and special printout statements for diagnosis and 
reminders. In each of the group descriptions below the first 
line contains the data label integer and the second the variable 
names . 


10 

DNUIC, DETAM 
where 

DNUIC - fractional allowable error on the Nusselt number 
calculation in Subroutine MARENU, Default value 
= 0 . 002 . 

DETAM - fractional allowable error on the recovery factor 

calculation in Program HOTWIRE. Default value = 0.002. 

11 

EBLT - an estimation of the boundary layer thickness. Default 

value = 0.10. This variable has no effect on any 
computations, all it does is control the distance over 
which a certain number (KK, see below) of data points 
are read; printout is "estimated boundary layer thick- 
ness = EBLT." 
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12 

DETAIL - a parameter to be set to non-zero if a detailed 

printout is required for diagnosis. Default 
value = 0.0 for no detailed printout. 

14 

NPROFL - the number of profiles or data sets to be reduced. 

used only if the flow direction is to be computed. 

No default value. 

15 

ITMARN , ITETAM 
where 

ITMARN - the maximum number of iterations allowed for sub- 
routine MARENU to determine the Nusselt number 
from the correlation curves. Default number = 10. 

ITETAM - the maximum number of iterations allowed for program 
HOTWIR to compute the recovery ratio. Default 
value = 5. 

Note that there are two controls on the required accuracy of Mach 
number and Nusselt number. If the maximum numbers of iterations 
are exceeded a printout "Solution not accurate" will occur along 
with the inaccurate solution. 

Tunnel and Wire Data 

These data specify the fluid, wire, and anemometer constants. The 
default values are for the fluid being air, and the wire a 0.0004 
inch platinum - 10% rhodium wire. The bridge resistance of the 
hotwire anemometer is 50 ohms. 
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GAM, R, PR, CMU 
where 

GAM - specific heat ratio. Default value * 1.4. 

R - gas constant. Default value * 53.61, 

PR - Prandtl's number. Default value = 0.71. 

CMU - fluid viscosity. No default. Computed internally 
using Sutherland's viscosity formula for air. If a 
different fluid is used, the Sutherland equation 
must be bypassed. 

5 

WD - wire diameter. Default value = 0.0004 inches. 

7 

RESIST - hotwire anemometer bridge resistance. Default 

= 50 ohms. 

8 

CALCON - parameter near 1 to calibrate wire diameter if 

necessary. See computational procedure step 25. 

Default value « 1.0. 

9 

WK, TK 
where 

WK - unit length resistivity of the wire at the given 
temperature. Default value = 57.25 ohms per inch. 

TK - temperature at which the above resistivity is 
specified. Default = 527. 7°R. 
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16 

WIKI, TW1 , WIK2, TW2 
where 

WIKI, WIK2 - thermal conductivity of the hot wire at two 

different temperatures. Default values = 43.35, 
42.47 BTU/hr-ft-°R. 

TW1 , TW2 - the corresponding temperatures. Default 
= 450, 720 °R. 

17 

ALFA, BETA 
where 

ALFA - temperature coefficient of resistance. Default 
= 0.96E-3 °R _1 . 

BETA - second temperature coefficient of resistance. 

Default = - 0.462E-7 °R -2 . 

18 

TT - total temperature. Default 535 °R. 

Specific Test Data 

These data include output voltages of the anemometer, temperature- 
resistivity calibration, static pressure, total temperature, overheat 
and wire lead resistances, and the wire orientation angles. 

1 K 

{Alpha-numeric title or blank) 

VOLTS(l), VOLTS (2) , V0LT$(3), ... , ... , V0LTS(6) 

V0LTS{7) , ... , VOLTS(I) , .... VOLTS(K) 
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where 

K - number of output voltages to be read for this 
particular data run, K < 200. 

(title) - any alpha-numeric or blank statement to be printed 
as a title for this test run. If omitted, program 
terminates. 

VOLTS(I) - the bridge voltage of the anemometer. Default a 0.0 
and program termination. 

2 

PINF - static pressure. Must be obtained separately as it is not 

part of the hot wire anemometry instrumentation. Default 
value = 0.495 psi . 

3 

Tl, WR1 , T2, WR2 
where 

WR1 , WR2 - the wire resistances as measured by the anemometer 
bridge at the corresponding temperatures. Default 
values = 2.6, 3.047 ohms. 

Tl , T2 - the corresponding temperatures. Default values = 492, 
672 °R. 

6 

WRWM, RESISL 
where 

WRWM - overheat resistance. No default. Omission causes 
program termination. 

RESISL - hot wire lead resistance as measured by the hot wire 
anemometer bridge . No default. 
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13 

THETAI , THETA1 , THETA2 
where 

THETAI - parameter to be set to 1 or 2 if the program is to 
determine local flow direction. If THETAI = l.,then 
data following Is for THETAI, and if THETAI = 2., data 
following for THETA2. Default = 0.0. If omitted, no 
angular calculation. 

THETAI - wire orientation or angle for the first Mach number 
component as shown in Figure 5. Default = 90°. 

THETA2 - wire angle for the second Mach number. Default = 45°. 

The execution starts after reading in VOLTS(K). Hence the data sets 
2, 3, 6, 13 must precede set 1. If THETAI is non-zero then two sets 
of data must be read in each with the same K. PINF need not be same 
but should not differ by more than 10%. The two data sets for the 
angles need not be together. A maximum number of profiles or data 
sets can be read in at one angle and then the profiles at the second 
angle may be read in. The only restriction is the number of profiles 
at each angle must be the same and the-sequence of profiles must be 
in the same order for each angle. 


7 
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5.0 CONCLUSIONS 

The data reduction procedure Is limited by the following restrictions: 

1) fluid can be gases only, which approximate the ideal gas 1aw» 

2) overheat ratio <.0.5, 

3) Reynolds number based on wire diameter 0.1 iRe^, <.300, and 

4) wire diameter should be calibrated. 

Other than these restrictions, the program is quite flexible. There are no 
restrictions on the Mach or Prandtl number. 

The default values in the computer program for the fluid and wire 
properties are those for air and a 0.0004 inch diameter platlnum-10% rhodium 
wire. These values may be changed if other gases and wires are used. The 
data input requirements are given in the previous section. The program 
listing Is in Appendix A. 
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FIGURE la. Schematic of Hot Wire in the Boundary Layer of the Model 



FIGURE lb. Coordinate System and Heat Losses of the Hot Wire. 
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FIGURE 2 


( b ) 


End loss correction for the Nusselt number ( a ) and 
recovery temperature ( b ), Reference 10. 
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FIGURE 3. Empirical convective heat transfer correlation for cylinders 
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EEGUKE 4 . Free-Stream Knuds en Number Correlation Curve, ( reference 10 ) 
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FIGURE 5. Coordinate System for the Local Flow Direction 
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APPENDIX A 

PROGRAM LISTINGS 


PROGRAM HOT'.') IR ( INPUT 4 OUTPUT 4 TAPE50 = INPUT) 
r 

c PPf'GP am hot WIPE the* computes mach nijmreps ano velocities 

c 

LOGICAL I EOF 

DIMENSION XMTH1 ( 20 ,prr>) , XMTH2 < 20 . ?00 ) 4 PINFI <?0> «PINF?(20) 
COMMON/RFS'JLT /XM 4 REI *CNUI »CKNl »cnum«cnuic 
COMMON/ETA/ ETA «ETAC »FTAF#ETAS » ETABAR ♦ FTAM 4 ET AM 1 ♦FTADEL 
COMMON/CCRCOE/ CN»CR 

COMMON/ENDLOS/ S * AWBAP 4 Al'/STAR 4 AWRA 4 OMEGA 4 OMGR A 
COMMON/TUNNEL/P INF »TT «TINF 

COMMON/") I PC AL/TV 1 4 W I !< 1 « T»i? « W I K P « T 1 «WR 1 « T? « WP? * WK' 4 TK I C ALCON 
COMMON /W I PF/ WD » WL 4 ALFA ♦FFT A « W I K I WREM 4 V'RWM 4 WRM 4 TFM * TWM 4 DFLT 
COMMON/GAS/GAM 4 PfPP 4 CMU,POAM 
COMMON/VOLTS/VOLTS < 200 ) 

COMMON /RES I ST/PFS I ST 4 PFS I SL 

COMMON/KONST/K 

COMMON/ I NOFX/ I J I 

COMMON/ INDEKS/KJI 4 JJJ 4 LLL 4 I TMARN 4 ITETAM 
COMMON/PREC I S/ONU I C 4 DETAM 
COMMON /EBLT/EBLT 
COMMON/DETA I L/OETA I L 

COMMON/THETA I /THETA I 4 THETA 1 4 THETA? 

COMMON /NPROFL/NPPOFL 

DATA ( VOLTS { I ) « 1=1 4?00)/?00*0.c/ 

DATA ALFA 4PFTA/.Q6F-3 4-4 A6PF-7/4 TK4('/K/F?7.7 4E7.?E/ 

DATA W IK1 ,TW1/43.3R44SO./,WIK24TW?/4?.4747?0./ 

DATA R,GAm,PR /53.61 * 1 *4i0.710/ 

DATA T I 4 WR1 4 T2 4 WR2 /49 1 • 7 1 2 • 6 4 67 1 • 7 4 3 . 047/ 

DATA PINE4TINF4TT/.495 4450.4535./ 

DATA WD 4 WL/O » 00C4 4 .0446/ 

DATA RESIST/50./ 

DATA RES I SL/. 24/ 

DATA DNU IC.OETAM/.OO? , .002/ 

DATA C ALCON/ 1./ 

DATA ERLT/ 0 . 10 / 

DATA DETAIL/0.0/ 

DATA (PINE1 < I) 4 1=1 4?0>/?O *0.0/ 

DATA (PINF2C I > 4 1=1 420)/P0 #0*0/ 

DATA THETAI 4 THETAI 4 THETA?/ 0.0490.445./ 

DATA ITMARN. ITFTAM/IO 45 / 


REPRODUCIBILITY OF THE ORIGINAL PAGE IS POOR, 


non non non 


1 «=J 


NPTHl s o 
NPTH? = 0 

CONTJMUr 

CAUL DATA{ IEOF) 

IF( I EOF > GO TO 5000 
IF ( TW=TAI .EQ. 1.1 GO 
IF ( THFTAI « FO » 2.) GO 
GO TO 4000 
CONTINUE 

N°TH 1 - NPTHl + I 


TO 

TO 


15 

16 


P INF1 ( NPTHl ) = PINF 

GO TO 4000 
16 CONTINUE 

NPTH? = NPTH?. + 1 
P INF? (NPTH?) = PINE 
4000 CONTINUE 

PGAM = (GAM - l.J/2. 
SGAMP I = SORT (GAM* I .5708) 
'--AMPO s t?. I 74*P*OAM 


OFTFPMINE A IP 
VISCOSITY ( 


THFQMAL GONOUC T I V/ 1 TY < RTU/HP-FT-R > 
LRM/FT-SEC) 


AND 


T 

= TT 


TIE 

= T**1.5 

1 

A KM 

= .001 1 40B*T 1 F>/ ( T + 

441 .72/ ( 10.*#(?1 .6/T) ) ) 

CMU 

= ( 7. 323E—7 ) *T 1 5/ ( T 

+ 198.72) 

CALIBRATE 


P 1 

= nip i 


R2 

= HR? 


ALFA 

= (R2-R1 >/(Rl*(T2-Tl 

> > 

t»'R 

= PI *< 1 .4-ALFA* (TK-Tl 

> > 

WL 

= I'.'R/WK 


PRINT 

?noo 


DETERM 

INF MFASURED HEATED W 

IRE TEMPERATURE 

TWM 

= T 1 + { < WRWM-WR1 )/(ALFA*WRl })-¥■(! • —RET A* ( WRWM— WR 1 

( WR1*ALFA#ALFA 1 ) 


WIK 

= W I K1 + ( TWM—TW 1 >*(WIK1-WIK2>/(TW1-TW2) 

KM 

= .OF 


ASSUME 

VALUE FOR THE ETAM 


ETAM 

* 1 . 



n r> 
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DO 1 O DO I a 1 » K 
JJJ = C 
7F TFM = PTAM-K-TT 
JJJ a JJJ + 1 

DEL.T = TWM-TCM 

WRFM = R 1 # ( 1 * + ( TEM-T 1 ) # ( ALF A + BFT A* < TFM-T 1 )) ) 

I'fRM = t,tp|nM 

F'-'P = WRFM 

VOLT = VOLTS ( n*WP , <iF/0' t RWM+RFSI ^T+RFSISL) 

OM = <VOLT*VOLT/’WRl>'M) 

CNUM = 13.04? * OM/( akm*WL*DELT > 

S = CV , D/WL>*SQRT<". | lK*WRM/<AI''M#FWR*CNUMn*CALCON 

awbar = ( WRM-EWR) /ewr 

AWRA = 1. + S*C-1 • 1+S*< .C999+S*<-?.?5) ) > 

?FO AWSTAR = AWRAP/AWPA 

ON = A"'Pfl#( 1 .+AWSTAP) /( 1 * + AUlRAP) 

CNUI = CN#CNUM 1 

I J I = I 

r 

r sol vf for mach numffr via itfration 

r 

nr j r = n 

call marfmu 

IF ( DETAIL • EO . O.O) GO TO 13 
IF ( I J I .NF. 1 .AND. IJI .NE. K> GO TO 13 
C #****#•#•*-•!* ■JHf'rt *****###****#*#*#*#*•** *■* ******* ******** 

PR I NT 10 

10 FORMAT (* ALFA.WL*WIK.WRM«EWR.TT«TEM /DELT ♦ ETAM . VOLT ♦ OM , CNUM . S . 
j AWBAR/AWRA. AWSTAR<CN»CNUI «CNUIC*REI < XM *) 

PRINT 1 2 « ALFA . WL ♦ W I K « WRM i EWR » TT » TF.M 
PRINT 1?. DFLT . FT AM * VOLT ♦ OM* CNUM ♦ S » AWBAR 
PPINT IP, AMIRA ♦ AI»iBTAR.rN.rNU I .ONUlCtRFI .XM 
1? FORMAT { 7F 14.4) 

C *##4f#*#*##*3Hf*#*tt*#R***#^#**##****###*******##**#*#**### 

1.3 CONTINUE 

CKNI = ( XM/RF I ) *SGAMP I 

RM3B = XM**3.5 

ETAC = 1 .-.05*RM35/( 1 . 175+RM35) 

RM28 = XM#*?.fl 

ETADEL = .2167*RM28X( .8F21+RM28) 

RKll = CKN I ## 1 . 1 93 
ETABAR = RK1 l/( .493+RK1 1 > 

ETA = FTARAR*ETADEL+ETAC 
RM2GAM = PGAM*XM*XM 

COMPUTE SUPPORT RECOVFRV RATIO a FTAS = TS/TT 

r 

ETAS = < 1 ,+SQRT(PR >*RM?GAM >/( 1 .+RM2GAM) 

OMGRA = S*< .7+S*(+l • 125+S#1 1 .25 >> 

QOO CONTINUE 

OMEGA = OMGRA/I 1 .-fOMGRA > 



non n o n n n n t p o 
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SOLVE FOR FNO LOSS PFC.OVFPY FAT. TOR CORRECTION PATIO 
CR = ( 1, -OMEGA* (FTAS/FTAM) )/< 1 .-OMFGA) 

COMPUTE FOR FTAM 
FTAM1 = FTA/CR 

COMPARE COMPUTFO AND ASSUMED ETA 

IF CARS! (ETAM-ETAMI J/FTAM) »LE. OETAM) GO TO 050 
9?0 CONTINUE 

JFIJJJ.GT, ITPTAM) GO TO 051 
FTAM = FT AMI — ( r T A M 1 — F TAM} /? . 

GO TO 75 
051 CONTINUE 
PRINT 718 

71 B FORMAT (* SOLUTION NOT ACCURATE- MORE. ETAM ITERATION NEEDED * ) 

950 CONTINUE 

I F ( LLL »GT • I TMARN ) PRINT 717 

717 FORMAT <* SOLUTION NOT ACCURATE - MORE MARFNU ITERATION NECSY *> 

COMPUTE LOCAL FREE STPFAM TEMPERATURE * DENS I TY 4 MASS FLOW* AND VELOCITY 


TINE = TT/( 1 4+RM2GAM) 

RHO = 1A4.*PINF/(P*TINF) 

V = XM*SQRT<GAMRC*TINF) 

IF ( DETAIL .EO. 0.0) GO TO 14 
IF { I J I .NF, 1 .AND. IJI *NE . K ) GO TO 14 
C ***W********************** 

PRINT 7001 

700 1 FORMAT <* CN i CNU I ♦ CKN I . ETAC tETADEL * ETABAR * ETA/ETAS 4 OMGRA 4 OMEGA , 
1 CR * ET AM 1 *T I NF * P I NF/WD 4 RE I 4 R 4RH0PT15 « V * RM2GAM #) 

PRINT 5001, CN.CNUI .CKNl * FTAC * ET ADEL ♦ ETABAR , ETA 
PRINT 5001, ETAS , OMGRA * OMEGA t CR 4 ET AM « T INF «P INF 
PRINT 5001, WD.REI.R. RH0*T15 ,V,RM2GAM 

5001 FORMAT (7*14.4) 

PRINT 6002 

600? FORMAT ( * TWM * DFLT 4 W I K » T ♦ CNU I C ♦ RE I * XM *> 

PRINT 5C0 1 , TWM 4 DELT 4 W I K 4 T 4 CNU I C 4 RE I 4 XM 

C ************************** 


14 CONTINUE 
A I = I 

PRINT 980 , A I 
9R0 FORMAT ( 1X#F3. 
?000 FORMAT { 1H 4 1 


1 

, 

1 OX , 2 OH 

?RE, 

R 

) 

IF 

( 

THETA I 

IF 

( 

THETA I 

GO 

TO 1020 


4 XM 4 V 4 RHO 4 T I NF 

0,F?0.4, 10x*F20.4 ♦10X4E20.4*10X*F?0.4) 

OX* 19H LOCAL MACH NUMBER 4 IOX 4 POH LOCAL VFLOC I TY 4 F/S 
LOCAL DENS I TY #L/CF 4 3X « 32 H LOCAL FREE STREAM TEMPERATU 

EO. 1 . ) GO TO 1015 
EQ. 2. > GO TO 1016 


REPRODUCIBILITY OF THE ORIGINAL PAGE IS POOR, 
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1015 CONTINUE 

XMTH1 (NPTH1 » I) = XM 
GO TO 1030 
10 18 COMTINUF 

XMTH?( NPTHP* I) = XM 
10p0 CONTINUE 

IF( (NPTHl+NPTH?) .FQ. <?*NPROFL) *AND. I.FG.K) GO TO 3000 
1111 CONTINUF 
LOGO CONTINUE 
GO TO 1 
3000 CONTINUE 

IF ( THETA I .NF. 1. .ANPiTHETAI .NF. 2«> GO TO 5000 
PTET1 = THETA1/57.2P6 
PTFT? = THETA2/57.296 
DO 8100 NP = ltNPPOFL 

IF ( PINFl(NP) .NET • P I NFH < NP ) ) GO TO 1017 
PINF = PINF1 (NP) 

GO TO 1018 

1017 CONTINUF 

PINF = (PINFl(NP) + PINFP(NP) >/?. 

1018 CONTINUE 
ANP = NP 
PRINT 8000 « ANP 

8000 FORMAT (1H1.F3.0,* PROFILF * > 

PRINT 8001 

8001 FORMAT ( 8X ♦ 32H LOCAL MACH NO. ANGLE ♦ DEG «24H LOCAL VELOCITY* 

1 F/3 » pX.PlH LOCAL DFNSITY* L/CF *2X*35H LOCAL FREE STREAM TEMP 

PERATURF * R ) 

DO PPOO J = 1 ,K 
AJ = J 

RTETA = ATAN( (XMTH1 (NP. J > #5 1 N ( RTET2 > - XMTH2 { NP ♦ J > *S I N ( RTET 1 ) > / 

1 ( XMTH 1 ( NP . J > *COS ( RTFTP > - XMTH2 C N° * J ) *C0S ( RTFT 1 ) ) ) 

THETA = PTETA*F7.PR6 

XM = XMTH I ( NP 4 J > /COS ( 1 • 5708 - RTFT 1 + RTFT A ) 

PMPfiAM = PGAM*XM*XM 

TINF = TT / ( 1 • + RM2GAM) 

RHO = 14 4* *P I NF/ ( R*T I NF ) 

V = XM*SQRT(GAMRC*TINF) 

PRINT 8002* AJ *XM* THETA* V« RHO «TINF 

8002 FORMAT < 1X»F3.P«E15.4»1?X»F5.1 . 3X . E 1 8 • 4 ♦ 7X ♦ EPO . 4 . 10X.EP0.4) 

8200 CONTINUE 

8100 CONTINUE 

NPTH1 = O 
NPTH2 = 0 
THETA 1 = 0.0 

THETA2 = 0.0 
GO TO 1 
5000 CONTINUE. 

STOP 

END 
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«*! ifpout lNr mapFMU 


CO^OM/Pr^ULT /XM 4 PFI 4 CMUI 40KNI ♦CMU^«CNU!C 

Common /FT A/ ft A »FTflP »P*TAr,rTAS 4 FTAPAP 4 FTAM 4 FTAMI »FTADFL 

COMMON/COPCOR/ CMtCR 

COMMPN/ENOLOE/ c * AWRAP 4 AV1ETAR 4 AWRA 4 OMFGA 4 OMGRA 
COMMON /TUMNFL/P ! NF 4 TT 4 T 1 NF 

COMMON /Vi I PCAL/TW 1 * V> l K 1 . TW2 4 W I K2 « T J 4 WR 1 4 T2 4 WR2 4 WK t TK 4 C ALCOM 
COMMON /'A I PF/WDf WL.4 ALFA 4 PET A 4 W1 K 4 WREM 4 WRWM 4 WRM 4 TFM 4 TWMtOFLT 
C OMMON/GA E/GAM 4 R 4 PR 4 CMU 4 RG AM 
COMMON/VOLTE/VOLTE < 200 ) 

common/komst/k 

COMMON/ I MDFX/ I J I 

COMMON/REE I FT/FFS I FT 4 RFE j cl 

COMMON/ I NDEKE/K J I 4 JJJ 4 LLL 4 ITMAPN 4 I TFT AM 

COMMON/PPF'C t F/ONU I C 4 OFTAM 

COMMON/FRLT/FRLT 

LLL = O 

YO = ALOGlO(CNUl) 

DLX = 0,1 

XI = ALOGIO(XM) 

KKK = 1 

GO TO 100 

1 CONTINUE 1 

I F( ARE < ( CNU IC-CNUI ) /CNUI C. 5 • LE • DNUfC.) GO TO 600 

Y 1 = ALOGlO(CNUIC) 

IOOO CONTINUE 

X2 = XI + OLX 

KKK = 2 

XM = 10,**X2 

CO TO 1O0 

2 CONTINUE 

LLL = LLL + 1 

I F ( LLL 4 GT * ITMARN ) GO TO 60 0 
C COMPARF CNU I C AMD CNUI AMD RETURN IF GOOD 

IFCAREC (CNUIC-CNUI 1/CNUIC) • LE 4 DNUICJ GO TO 600 
Y2 = ALOG 1 0 ( CNU I C ) 

DLX = (YO - Y2)*(X2 - XI >/(Y2 -Y1 ) 

XI = V? 

Y 1 = Y2 

GO TO 1000 
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ino CONTINUE 

PEI = 6a.06655*(XM#P!NF*WD/CMU)*CALC0N* 

1 SORT ( GAM# ( 1 • + RGAM*XM#XM)/(R*TT ) ) 

C IF XM LFSS THAN .4 THEN AN AVERAGE IS TAKEN OF THE VALUES OF CNUIC 
r at XM = 0.0 aMO XM = ,4 

IPIXW-.4) 400 » 400 , 4 ?jO 
400 CONTINUF 

Dfi = PFt**0,6713 

5N = 1 i-(R6)/(F. 14?+?**R6) 

R7 = RE I **0.71 14 

CN1 = 0.?30?*R7/< 15.44+R7) 

CNF = ( .0l569/< .3077+REI**.737B>>*< 15./I 1F5.+RFI**3. > > 

PH II = 1.00283 

Rll = RE I **1.109 

PHI? = 1.834 - 1 ,634*{R] 1/I2.765+RJ 1) ) 

PH 1 3 «1.- .OO0?42*(REI/(4. + RED) 

PHI = 1. + PHI 1*PHI?*PH13 

CNUIM4 = ( . 14 + 0.M1+CN? )*(RFI**SMl*PHI 

CNUIMO = .19 + .66*PFI**.45 

FF = -.618034 + 1 ./( (XM/.4 )**.9 + .618034) 

CNUIC = 10,**(FF*ALOG10(CNUIMO) + f 1 • - FF ) *ALOG 1 0 < CNU I M4 ) ) 

GO TO 500 

450 CONTINUE 

P6 = RF I **0 .67 1 3 

SN = 1 .-<R6>/<5. 142+2.*R6) 

R7 = RE I **0 .71 14 

CN1 = 0.2302*R7/( 15.44+R7) 

CN2 = ( .0156R/( .3077+RFI**»7378) )*( 15./( 15.+REI**3. ) ) 

PM 1 ? = XM** 1 .??? 

PH I 1 = .6039/XM + .5701*1 (RM12/U • + RM 1 2 ) > ** 1 .569 - 1.) 

Rll = RE I **l«109 

PHI? = 1.834 - 1 ,634*<RH/(2.76«+Rl 1) ) 

PH I 3 = 1 . + ( ,3-.065/<XM**l .67) > * < RE I / < 4 . +RE I ) > 

PHI = 1. + PHI 1*PHI?*PHI3 

CNUIC = C . 14+CNI+CN2 >*( RF I **SN > *PH I 
500 CONTINUE 

GO TO (1.2J.KKK 

600 RETURN 
FNn 
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SUBROUTINE DATA(IFOF) 


LOGICAL IFOF 

COMMON/RESULT /XM»REI »CNUI *CKNI *CNUM«CNUIC 

COMMON /ET A/ FTA * ETAC ,FTAF * ET AS * ET A BAR « FTAM, ETAM l ,FTADEL 

COMMON/COPCOF/ CNtCR 

COMMON/FMOL.OF/ S » A'i'RAR* A |,f< ^TAR» AWRA « OMFGA » OMGRA 
COMMON/TUMNFL/P I NF ♦ TT , T I NF 

COMMOM/WIPCAL/TWI ,WIK1 ,TWP* WIK?»T! ,WR1 * T? ♦ WR? * WK » TK , CALCON 
COMMON/'*' I PF/ WO ♦ WL , ALF A , BETA * W I K « WRFM « WRWM « WPM » TFM , TWM , DEL T 
COMMON/GAF/GAM < P « PP « CMU < PGAM 
COMMON/VOLTS/VOLTS (200 ) 

COMMON/KONST/K 
C0MM0N/RFS I FT/RES I ST ♦ RFS ! SL 
COMMON/INDEKS/KJI , JJJ.LLL* I TMARN * ITETAM 
C OMMON/PRFC I S /DNU 1 C ♦ DETAM 

common/eblt/eblt 

COMMON/DETA IL/DETA I L 

COMMON/THFTA ! /THFTA I » THFTA I * THFTA2 

commom/nppofl/nprofl 


1000 COMTINUF 

OBAn 1 00 , n »KK 
100 FORMAT ( IP * IS) 

IF ( EOF ♦ SO ) 2000* 1 01 
10 I CONTINUE 

I EOF = .FALSE. 

I F (KK.GT, 0) K = KK 
C 

C KK- SHOULD RE READ IN ONLY IF II s I AND IF NUMBER OF 
C VOLTAGE READINGS PER RUN CHANGES 
C 

GO TO ( 1 ,2, 3* A , B,6» 7.B.O* 1 0. 1 I « 1 2 « 1 3t 14* 1 F* 16. 17* 1 fl « I R» 20) » I I 
C 

c note all units arf as follows 
c resistances ohms 
c temperatures degrfes rankine 

c lfngths inches 
c viscosity lbm/ft-sec 
■ c gas constant ft-lbf/lbm-r 

*c conductivity btu/hr-ft-r 

c 

C A card must prfcfdf each data block to lapel the following data 

c. The card prfcfdjng the voltage data contains two inputs. II AND k 

C The label card for all other data blocks contains only one input, ii 

c 




1 
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c 

c 

c 

c 

c 

c 

c 

c 

r 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

C 


II a J SHOWS THE VOLTAGES OF ANFMOMFTFP READINGS 
K * NUMBER OF voltage pfadings TAKEN for THAT RUN 


I 1 * 

II *s 


I I = 


I I = 
I I = 

I I = 

I I = 

I I = 


P LOCAL EopF STREAM STATIC PRESSURE , PINF 
3 V/IRF CALIBRATIONS T I ♦ WR 1 t T2 , WR2 
A GAS DATA GAM*R*PR*CMU 
F WIRE DIAMF1FR* V'D 

6 HFATFD WIRE RFS I STANCE . WRWM * AND LEAD RFS I STANCE » RES I SL 

7 FRIDGE rf SI STANCE PS t RESIST 
B CAL I DRAT I ON CONSTANT* CALCON 

O WIRE UNIT LENGTH RFS I ST I V I T Y < OHMS/ 1 N * > «WK « AT TEMP.* TK 
P PRFCISrON ON NUSSFLTG NO. AND EO. TEMP. RAT I O * DNU I C .DFTAM 
J1 ESTIMATED n.L. THICKNESS* EBLT * (IN.) 

1? PARAMETER TO BE SFT TO NON-ZERO IF A DFTAILFD PRINTOUT IS 

REQUIRED, DFTAIL. THE P.O. WILL HE FOR THF FIRST AND LAST 
RUNS ONLY. 

13 PARAMETER TO BE SFT TO I. OR 2. IF ANGLE DEPENDENCE IS REOD 
THFTAI = 1. FOR 1ST ANGLE 

THETA I s ?. FOR 2ND ANGLE 

AND THE VALUES OF THOSE TWO ANGLES* THETA 1 « THFTA2 

14 NUMBER OF PROF ILFS TO RE READ IN AT EACH ANGLE, NPROFL 

15 NUMBER OE ITERATIONS TO COMPUTE NU*ITMAPN* AMD COMPUTED 
RECOVERY RATIO, ITFTAM 

16 THERMAL CONDUCTIVITIES OF WIRE AT TWO TEMPERATURES, 

WIKI *TW1 *WIK2»TW? 

17 FIRST AND SECOND TEMPERATURE COEFFICIENT OE RESISTIVITY, 

ALFA AND BETA 

IB TOTAL TEMPERATURE , TT 


ALL DATA MUST BE RIGHT JUSTIFIED 


1 DPAH 200 

200 format (70H 

1 ) 
IF ( EOF, RO ) 2000,201 

201 CONTINUE 

I EOF = *FALSF. 

PRINT 300 
300 FORMAT ( 1 HI ) 

PRINT 200 
PRINT 350 
350 FORMAT (1H > 

READ 401 , (VOLTSU) ♦ 1 = 1 *K> 

401 FORMAT ( FB, 4 * 5F 1 0 • 4 > 

400 FORMAT (6F10.4) 

RETURN 

2 READ 400, PINF 
GO TO 1000 

3 READ 4 CO , T 1 , WR 1 » T2 • WR2 
GO TO 1000 

4 READ 400, GAM , R ,PR « CMU 
GO TO loon 



REPRODUCIBILITY OF THE ORIGINAL PAGE l$ POOR, 
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5 READ 400 ♦ WD 
GO TO 1000 

6 REAR 400 « WPV/M»RESISL 
GO TO JOOo 

7 RFAD 400, RFSIPT 
GO TO POO 

B PEA O A 00, CALCON 
GO TO lOOO 

P READ ' 4 nr| , WK»TK 
GO TO 100o 

10 REAO 400 , DMU I C »DFTAM 
on to iooo 

11 READ 400, FPLT 
PRINT 9B2 , F«LT 

OB P FORMAT ( * FSTIMATFD BOUNDARY LAYER THICKNESS = *»F6.3« * INCHES* ) 

GO TO 1000 

1 P READ 40-, OFTAIL 
GO TO 100n 

13 READ 400, THFTAI ,THFTAI 4THETA2 
GO TO 1000 

14 PEAO ICO, NPROFL 
GO TO 10C0 

15 READ 100, 1 TMARN , I TFT AM 
GO TO 100Q 

16 READ 400, V IK1 »TW1 »W IKP,TW2 
GO TO 1000 

17 READ 400, ALFA,RETA 
GO TO 1^00 

IB READ 400 , tt 
GO TO 1000 

IQ CONTINUE 
GO TO 1000 

20 CONTINUE 
GO TO 1000 
2000 CONTINUE 

I EOF = .TRUE, 

RETURN 

END 
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APPENDIX B 

NON-LINEAR RESISTIVITY EFFECTS 


The effects of neglecting the second temperature coefficient of 
resistivity are considered in this appendix. 

The resistance of the wire at temperature T w is 

r w ■ r oO + “l<V T o> + «2<W Z] 

where 

r Q = resistance at temperature T Q 

a 1 9 first temperature resistivity coefficient (linear) 

a 2 = second temperature resistivity coefficient. 

Rewriting in terns of overheat ratio 


r -r 
w o _ 


* “1<V T 0> ww 


and solving for (T w -T 0 ) 


t t - J_ r w~ r o 

w " 0 ’ a l >0 


if cu = 0 


r r - 1 ®1 m * 4a 2 V p o J/2 n a n 

T w ~ T o 7 57 W + 5T ”r“ ) “ 13 lf a 2 * 0 

2 1 0 

expanding the bracketed term with 


i 4o 2 , Vo 
“T ( -r~~ 

Ot-j 0 


< 1 


yields 



50 


r t - 1 Vr ° 1 ^ f 

w o ' “! r 0 " T r 0 


Qo 


w" r o %Z . 1 a 2 , V r o x3 


) + 


37 ^ 


( 


r + ... 


therefore the error in the temperature due to neglecting the non-linear 
effect a 2 is 


a 1 T = 


r -r 
w o 


1 

T 


° 2 (^) 


1 V , [rt ,2 + 

37 ~T 1 ^ * 


a 


1 


-7 

The temperature error for typical values of |a 2 | 10“ and ^ ^ 10“ 

(Ref. 2) is 1% for an overheat ratio of 0.4 and 2.5% for an overheat 
of 1.0. 

The non-linear i&mperature coefficient of resistivity may be 
safely neglected if the overheat ratio does not exceed 0.4. 
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APPENDIX C 

UNCERTAINTY ANALYSIS OF THE MIRE DIAMETER 

The effects of the wire diameter uncertainties are considered 
here. The uncertainties of all the required input data should be con- 
sidered in determining the uncertainty of the output of the computational 
procedure of Section 3.0, i.e., Mach number. However, the diameter of the 
wire is the most uncertain of all the data inputs, and furthermore, the 
diameter is not readily measurable, either directly or Indirectly. Usually 
the manufacturer's nominal diameter is used {accuracy ±3%). The remaining 
data can be determined rather accurately, hence their uncertainties will 
not be considered. 

The uncertainty of a function, Af/f (e.g., Mach number), is 
defined by Kline and McClintock (13) as 


Af 

T 


- c l ff ) 2 (^) 2 ] 1/2 

i=l f 3x i x i 


= C l (x 1 

i=l 1 


9 an f x2 


9X 


1 


AX- 

% 


)‘ ( — ) 2 ] 1/2 


where f = f(x^), and the x.-'s and (Ax^/x^J's are the independent data 
inputs (and their respective uncertainties) required to determine f. 

The uncertainty of mass flow due to uncertainty of wire diameter 
is determined as follows 


Nu = A (Re) 1/2 + B 


where for 

incompressible flow A = 0.690 B = 0.318 

(King's Eq., Ref. 1) 
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and 

supersonic flow 


A = 0.580 


B = -C 


Nu = 


_ hd 


Re = 

n 


From the defining uncertainty equation 


_ / hd B i2jj 

e u - ( W - 5 J S' 


£n pU = JLn p - Jin d + 2£n ( £ ” ^ ) 




A 8 An pL) hd + Bk _ Nu + B 
Q " 9cf " ” hd - 6k ” Nu "HT 


*8 ■ t( % ) z ♦ ( £41 ) 2 ( “ ) 2 f /2 • 


For typical values in the subsonic range 

( M ) = ± [( 0.01 ) 2 + ( | - H {| ) 2 ( 0 . 10) 2 ] 1/2 

pu M<<1 ‘ - U.JIO 

~ ± 14% if = ± 10% 

and M = ± i% 

r* 

and Nu = 2 . 

For typical values in the supersonic range 


y 


.795 
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* * t(°.°i ) 2 + ( jprogf ) 2 (o-io ) 2 

~ ±8% If ^ s ±10% 

and ^ = ±1% 

and Nu = 8 

Therefore the uncertainty in the mass flow is of the same order as the 
uncertainty in the wire diameter. Since the mass flow is directly 
proportional to the Mach number, the uncertainties in Mach number have 
the same values. Hence an effort must be made to correct for uncertain 
wire diameters. 


